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Abstract. We discuss the creation of entanglement between two two-level atoms in the dissipative 
process of spontaneous emission. It is shown that spontaneous emission can lead to a transient 
entanglement between the atoms even if the atoms were prepared initially in an unentangled state. The 
amount of entanglement created in the system is quantified by using two different measures: concurrence 
and negativity. We find analytical formulas for the evolution of concurrence and negativity in the system. 
We also find the analytical relation between the two measures of entanglement. The system consists of 
two two-level atoms which are separated by an arbitrary distance ri2 and interact with each other via 
the dipole-dipole interaction, and the antisymmetric state of the system is included throughout, even for 
small inter-atomic separations, in contrast to the small sample model. It is shown that for sufficiently 
large values of the dipole-dipole interaction initially the entanglement exhibits oscillatory behaviour with 
considerable entanglement in the peaks. For longer times the amount of entanglement is directly related 
to the population of the slowly decaying antisymmetric state. 
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1. Introduction 

Entanglement is a property of quantum systems to exhibit correlations that cannot be accounted for 
classically. Entangled states of collective quantum systems, which are states that cannot be factorized 
into product states of the subsystems, are of fundamental interest in quantum mechanics. A number 
of methods for creating entanglement have been proposed involving trapped and cooled ions or neutral 
atoms m^lElBElEllZIIHl • Of particular interest is generation of entangled states in two-atom systems, 
since they can represent two qubits, the building blocks of the quantum gates that are essential to 
implement quantum protocols in quantum information processing. It has been shown that entangled 
states in a two-atom system can be created by a continuous driving of the atoms with a coherent or 
chaotic thermal field [511511 1U|. or by a pulse excitation followed by a continuous observation of radiative 
decay Moreover, the effect of spontaneous emission on initially prepared entangled state 

has also been discussed ^JEIEIE]- These studies, however, have been limited to the small sample 
(Dicke) model |18j or the situation involving noninteracting atoms strongly coupled to a cavity mode. The 
difficulty of the Dicke model is that it does not include the dipole-dipole interaction among the atoms and 
does not correspond to realistic experimental situations of atoms located (trapped) at different positions. 
In fact, the model corresponds to a very specific geometrical configuration of the atoms confined to a 
volume much smaller compared with the atomic resonant wavelength (the small-sample model). The 
present atom trapping and cooling techniques can trap two atoms at distances of order of a resonant 
wavelength jl9l2()l2i] . which makes questionable the applicability of the Dicke model to physical systems. 

Recently, we have shown [22] that spontaneous emission from two spatially separated atoms can lead 
to a transient entanglement of initially unentangled atoms. This result contrasts the with the Dicke model 
where spontaneous emission cannot produce entanglement from initially unentangled atoms jl0[ll6| . We 
have numerically calculated the evolution of the concurrence and discussed the role of the maximally 
entangled collective states of the two-atom system: the rapidly decaying symmetric state and the slowly 
decaying antisymmetric state of the two-atom system. 

In this paper we extend our study of spontaneously induced transient entanglement in a system 
of two atoms separated by an arbitrary distance ri2. We find analytical results for the two calculable 
measures of entanglement concurrence and negativity establishing the relation between the two. Our 
solutions are valid for a broad class of initial conditions including mixed states. It is shown that when 
the dipole-dipole interaction becomes larger than the atomic decay rate then the entanglement exhibits 
oscillatory behaviour, oscillating with twice the frequency describing the dipole-dipole interaction, which 
is the frequency separation between the symmetric and antisymmetric states of the two-atom system. 
Remarkable amounts of entanglement can be obtained at the maxima of the oscillations. For times longer 
that the decay rate of the superradiant symmetric state, when the population of the symmetric state is 
negligible, the oscillations disappear and the entanglement remaining in the system is related to the 
population of the slowly decaying antisymmetric state. 

2. Measures of entanglement 

To assess how much entanglement is stored in a given quantum system it is essential to have 
appropriate measures of entanglement. A number of measures have been proposed, which include 
entanglement of formation [53] , entanglement of distillation [21] , relative entropy of entanglement [^ and 
negativity [5B','77"3S','25" . For pure states, the Bell states represent maximally entangled states, but for 
mixed states represented by a density matrix there are some difficulties with ordering the states according 
to various entanglement measures; different entanglement measures can give different orderings of pairs 
of mixed states and there is a problem of the definition of the maximally entangled mixed state '30''31J . 

Here we use two entanglement measures, i.e., concurrence and negativity to describe the amount of 
entanglement created in a two-atom system during spontaneous emission. The concurrence introduced 
by Wootters [23j is defined as 




(1) 
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where {A^} are the the eigenvalues of the matrix 

R^pp (2) 

with p given by 

p = da; (g) CTa; p* O-^ (g) CTj; , (3) 

ax is the Pauh matrix, and p is the density matrix representing the quantum state. The range of 
concurrence is from to 1. For unentangled atoms C = whereas C ~ I for the maximaUy entangled 
atoms. 

Another measure of entanglement we use here is the negativity, which is based on the Peres- 
Horodecki |26[I77| criterion for entanglement and is defined by the formula 



M = max 0, 



(4) 



where the sum is taken over the negative eigenvalues p,i of the partial transposition of the density matrix p 
of the system. The partial transposition means transposition with respect to the one atom only. For pure 
states, like in the case of concurrence. A/" = 1 for maximally entangled state and JV = for unentangled 
atoms. 

The two entanglement measures, i.e., concurrence and negativity, give the same criteria for 
entanglement, but generally they give different values for a degree of entanglement [31] ■ We will give 
analytical expressions for both of them for the entanglement produced in spontaneous emission. 

Introducing the computational basis for the two-atom system as product states of the individual 
atoms, as follows 

|l) = l5i>'^l52) , 

|2) = |ei)® |e2) , 

|3) = \9i) ® |e2) , 

|4) = |ei) ® 132) , (5) 

where \gi) and je^) (for i — 1,2) are the ground and excited states of the individual atoms, we can define 
the density matrix of the two-atom system as a 4 x 4 matrix. 

We assume that the density matrix of the system has the block form 

/ Pll Pi2 ^ 



(6) 



P21 P22 
P33 P34 
\ P43 P44 / 

with the condition Trp = 1. We will show that, if the atoms initially start from a state described by 
the density matrix of the form the evolution does not destroy this form in the sense that the blocks 
of zeros remain untouched. The other matrix elements evolve in time, and we find explicitly their time 
dependence. 

The matrix p, needed for calculation of the concurrence, has the form 

f P22 Pl2 ^ 

P21 Pll 





P44 P34 
Pi3 P33 ) 

and the square roots of the eigenvalues of the matrix R given by ((2Jl are the following 



V 



(7) 



(8) 



|V^} = {\/PllP22 - \pi2\,y/pilP22 + |pi2 | , \/P33P44 - |p34 | , \/P33P44 + |p34| } . 

Depending on the particular values of the matrix elements there are two possibilities for the largest 
eigenvalue, either the second term or the fourth term in (jSJ. The concurrence is thus given by 



C = max{0, Ci, C2} 



(9) 
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with 



Ci = 2 (IP12I - VP33P44 ) , 

C2 = 2(|p34|-VpTTp^) ' (10) 
and we have two alternative expressions for the concurrence depending on which of them is positive. 

For calculation of the negativity we need the partially transposed density matrix. The transposition 
with respect to the indices of the first atom gives the matrix 



/ Pll 


P43 







P34 


P22 














P33 


P21 


V 





P12 


P44 / 



(11) 



which has the eigenvalues 

^ (p33 + P44 ± \/ {P33 + PiiY + 4 (|pi2p ^ P33P44)) 



Pll + P22 ± V (Pll + P22f + 4 (I/O34P - P11P22) 



(12) 



There are two candidates for being negative among the roots (|12|l . however, they cannot be negative 
simultaneously because the inequality |p34| — ^ p\\P22 > implies that |pi2| — ^/p33P4A < 0, and vice 
versa. It is also easy to find the two alternative values for the concurrence @ inside the square roots. 
So, the negativity defined by (@J) has also two alternative forms 



P12I 



/033P44j 



iP3 



Pii 



V4(|p34p - P11P22) + (Pll 
0, J Ci + (/733 + p44)^ - {P33 + P44 



{P33 + P44 
P22Y - (pil 



P22 



\IC2C2 + (/On + /922)^ - (/Oil + P22) 



(13) 



(14) 



where the appropriate expression from 1)101) is to be substituted to the two alternative terms in H14|l . The 
quantities Cj*", and C^, which are always nonnegative, represent two alternative expressions associated 
with corresponding expressions Ci and C2 for the concurrence lO, and they have the following form 

Cj^ = 2 (IP12I + VP33P44), 

C2+ = 2(|p34| + VpITp^)- (15) 
For pure states and are equal to Ci and C2, respectively, and in this case the negativity is equal 
to the concurrence. 

The equality H14|l establishes the relation between the negativity and the concurrence for the system 
described by the density matrix of the form It is evident that both quantities give the same criterion 
for entanglement, that is, positive value of C implies positive value of M, but the degree of entanglement 
indicated by the two quantities can be quite different. It is also clear from (fT^ that the product CiC^ 
(C2 C2) can as a whole serve as a measure of entanglement: it is zero if C — 0, and it is unity for maximally 
entangled pure state, but for mixed states it gives yet another value for the degree of entanglement. 

It is interesting to express the results for concurrence and negativity in the Bell basis which is defined 
as follows 



in 



12') 



i$+> 
i$-> 



13') = |*+) = 



-|2» 
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14') 



1 



V2 



(|3) - |4)) 



(16) 



The transformation of the density matrix p given by © from the original basis ^ to the BcU basis H16(l 
is performed with the transformation matrix 

/ 1 1 0^ 




1 



V2 



V 





1 1 

-1 1 / 



-1 1 




leading to the new density matrix 
p' = UpU+ , 

which has the same block form as (jSJ but the new matrix elements have primed indices {pvi 
The matrix elements in the Bell basis are related to the original matrix elements as follows 

Pl'l' = 5 [Pll + P22 + (/3l2 + P2l)] , P3'3' = 5 [/O33 + P44 + (Pl2 + P2l)] , 



(17) 



(18) 



,Pl'2' 



(19) 



P2'2' = 2 + '^22 - (Pl2 + P21)] , P4'4' = j [P33 + P44 - (/534 + P43)] 

Pl'2' = "2 [Pll - P22 + {Pl2 - P21)] , P3'4' = "2 [P33 " P44 + (P34 - P43)] 

P2'l' = -2 [Pll - P22 - (Pl2 - P21)] , P4'3' = -2 [P33 - P44 - (P34 - P43)] ■ 

In the Bell basis (|lt)|l . the concurrence alternatives (|10|l and the negativity H14|l take the following form 

Cl = V (Pl'l' ~ P2'2')^ ^ (Pv 

(20) 



P3'3' 



- P2'v f - V (/03'3' + P4'4')^ - (P3'4' + P4'3')^ 
■ P4'4')^ ' (/03'4' - P4'3')^ - V (Pl'l' + P2'2')^ - (Pl'2' + PTvY 

0, \JCiC^ + (/93'3' + Pi'A'Y - (P3'3' + P4'4') , 



N ■ 



C2C} + {pvy 



+ P2'2>Y - (Pl'l' + P2'2') 



(21) 



where 



P2'2' 



P2'vY + \J (P3'3' + PVvY - {p-i'V + P4'3')' 



V (/03'3' - PA'VY - (P3'4' - PV-i'Y + V (Pl'l' + P2'2')^ ^ (Pl'2' + /02'1')^ 



(22) 

From (|2()|l and (|21|1 . it is evident that for any Bell state H16|l the concurrence and negativity become unity. 
For mixed states the situation is much more complicated with the values of concurrence and negativity 
which are different in this case, and have values between zero and unity. Later on we apply the general 
formulas derived in this Section to find the evolution of entanglement in spontaneous emission from a 
system of two two-level atoms. 



3. Atomic evolution 



We consider a system of two non-overlapping two-level atoms with ground states \gi) and excited states 
\ei) {i = 1,2) connected by dipole transition moments fli. The atoms are located at fixed positions ri 
and f2 and coupled to all modes of the electromagnetic field, which we assume are in the vacuum state. 
We consider spontaneous emission from identical as well as non-identical atoms prepared in different 
initial states. In the case of nonidentical atoms, we assume that atoms have equal dipole moments 
Ml = ^2 = but different transition frequencies lui and lj2, such that 0J2 — oJi ujq = {oji + a;2)/2, so 
that the rotating-wave approximation can be applied to calculate the dynamics of the system. 

The time evolution of the system of atoms coupled through the vacuum field is given by the following 
master equation pT^r^lHl] 

~^j2^,[S!,p]-^j2^.AStS-,p] 
1 2 

- 2 E {pStS- + S+S-p - 2S-pSt) , (23) 



Entangling two atoms via spontaneous emission 



6 



where (S^) are the dipole raising (lowering) operators and is the energy operator of the ith atom. 
In Eq. H23|) . {i — j) are the spontaneous emission rates of the atoms, equal to the Einstein A coefficient 
for spontaneous emission, whereas Tij and flij {i ^ j) describe the interatomic coupling |32[l33ll5^ . and 
are the collective damping and the dipole-dipole interaction potential defined, respectively, by 

sin(fcory) 



r- — r ■ — 



f - (M • Ti, 



1 - 3 (/i • fj 



kortj 
cos(fcor^j) 



sin {kprij) 



(24) 



and 



1^ 



1 - (A ■ f^^jy 



f - 3 (/i • f-ijY 



COS {kprij) 
sin(fcorij) cos (kar^j) 



(fcory) 



(25) 



where ko — loq/c, Vij = \rj — ri\ is the distance between the atoms, p, is unit vector along the atomic 
transition dipole moments, that we assume are parallel to each other, and fij is the unit vector along the 
interatomic axis. 

The master equation has been used for many years to study a wide variety of problems involving 
the interaction of collective atomic systems with the radiation field [321 ■ Using the master equation (|23|l , 
we can write down the equations of motion for the components of the density matrix of the two-atom 
system in the basis (|5)l of the product states of the individual atoms. However, the problem simplifies by 




Figure 1. Collective states of two identical atoms 



working in the basis of the collective states of the system which contains symmetric and antisymmetric 
combinations of the product states. For identical atoms (wi = ^^2) the collective states are [THIIS^ 

\g) = |i) , 

|e) = |2) , 

k> =^(|3) + |4)) , 

|a) = i=(|4)-|3)) , (26) 
where we used the basis (jSJ. 
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In the collective state representation, the two-atom system behaves as a single four-level system, 
illustrated in Fig. ^ with the ground state \g), the upper state |e), and two intermediate states: the 
symmetric \s) and antisymmetric \a) states. The most important property of the collective states is 
that the symmetric and antisymmetric states are maximally entangled states. The states are linear 
superpositions of the product states which cannot be separated into product states of the individual 
atoms. They are in fact two of the Bell states introduced in H16|l : |s) = |3') and \a) = |4'). The 
symmetric and antisymmetric states are eigenstates of the system of two identical atoms with the dipole- 
dipole interaction included. The basis of atomic states (|26|l can be considered as an effect of partial 
transformation to the Bell basis IjKil) in which the transformation has been performed in the lower block 
only. Such basis is convenient for finding the solution to the master equation H23I) describing spontaneous 
emission in the system. 

Assuming that initially the state of the system has been prepared in the block form jnj, from the 
master equation (|^ . we get the following set of differential equations describing the evolution of the 
system in the basis of collective atomic states (|26f) 

Pee — 2rpee , 

Peg = - (r + 2iuJo) Peg , 

Pss = - (r + ri2) {pss - Pee) + «A {pas ^ Psa) , 

Paa = - (r - ri2) {paa - Pee) ~ «A {pas - Psa) , 

Pas = - (r + 2if7i2) Pas + «A (p^s - Paa) (27) 

with the condition pgg + pee + Pss + Paa = 1, and with A = (a;2 — wi)/2. All other matrix elements related 
to the blocks of zeros in ^ remain zeros, if the evolution is govern by the master equation 1)23(1 . 

Equations ((27|l show that all transitions rates to and from the symmetric state are equal to (r + ri2). 
On the other hand, all transitions rates to and from the antisymmetric state are equal to (F — ri2). Thus, 
the symmetric state decays with an enhanced (superradiant) rate, while the antisymmetric state decays 
with a reduced (subradiant) state. Hence, the population of the antisymmetric state experiences a 
variation on a time scale of order (F — Fi2)^^, which can lead to interesting effects not observed in the 
Dicke model. These effects result from the fact that the set of equations ((27|l has two different solutions 
depending on whether F12 = F or F12 7^ F. The case of F12 = F corresponds to the small sample (Dicke) 
model, whereas the case of F12 7^ F corresponds to spatially extended atomic systems. The existence 
of two different solutions of Eq. (|27|) is connected with conservation of the total spin S'^, that S'^ is a 
constant of motion for the Dicke model and S'^ not being a constant of motion for a spatially extended 
system of atoms [HEIEZI- We can explain it by expressing the square of the total spin of the two- atom 
system in terms of the density matrix elements of the collective system as 

{t) = 2~- 2paa (t) . (28) 

It is clear from Eq. ((28|l that S'^ is conserved only in the Dicke model, in which the antisymmetric state 
is ignored. For a spatially extended system the antisymmetric state participates fully in the dynamics 
and S'^ is not conserved. The Dicke model evolves between the triplet states |e), |s), and \g), while the 
spatially extended two-atom system evolves between the triplet and the antisymmetric states. 

The problem of spontaneous emission from two atoms can be solved analytycally even for general 
case of nonidentical atoms (A ^ 0) |SH1) but the general solutions are rather lengthy, and we will give 
here the solutions for the simpler case of identical atoms only. It is seen from 1(27(1 that the first two 
equations, belonging to the upper block of ©, are decoupled from the other equations belonging to the 
lower block of ©, and they have simple exponential solutions 

Pee(i) =Pee(0)e-"" , 

Pe<,(t)-Peg(0)e-(i^+2-")* . (29) 

For identical atoms, A = 0, the remaining equations simplify considerably and their solutions are as 
follows 

p^sit) = PssiO) e-(r+r..)* + ^^^(0) £±£ii (e~(r+r..)* _ ^-2Tt\ 

i — i 12 ^ / 
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i + i 12 ^ / 

The evolution within the two blocks runs independently except for the fact that all the states decay to 
the ground state \g), and the population of this state is 

Pggit) = 1 - Pee{t) - Pssit) - Paa{t) , (31) 

i.e, eventually total atomic population accumulates in the ground state. 

The solutions H29I) and H3U|I are particularly simple in the basis (|26ll . but it is quite easy to 
transform them into the original basis Q or into the Bell basis (I16II using the relations (|19() . The 
same transformation can be used to transform the initial conditions. In this way we obtain analytical 
results for the matrix elements of the density matrix in either the original basis or the Bell basis as 
the linear combinations of the solutions H29|l and (|30|l . for any initial conditions that preserve the block 
form © of the density matrix. 



4. Entanglement in the two-atom system 

The solutions obtained in the previous Section can be used in formulas l|10|l or H2Q|I for the concurrence 
and in formulas (|13|l or H21|l for the negativity giving the analytical expressions for the quantities 
describing degree of entanglement in the system. For example, if Peg(O) = P2i(0) = 0, i.e., there is 
no two-photon coherence in the system initially, then Ci cannot be positive, so it cannot contribute to 
the concurrence C, and the concurrence is equal to C2, if it is positive. We have 

C{t) = max{0,C2(t)} , (32) 

and 

C2{t) = ^[Pssit) - Paait)f " [psa(t) " Pasit)f ~ 2^ pee{t)pgg{t) (33) 

with the solutions (|23l-lEIJ- The solution (O stiU covers a broad range of initial conditions, i.e., such 
that the upper block in © is diagonal but the lower block is arbitrary. 

It is immediately seen from (I33|l and the solutions (|3U|I that, for two identical atoms prepared initially 
in one of the maximally entangled states \s) or |a), the concurrence for any time is equal to the population 
of the corresponding state Pss(t) or Paa{t). it is unity at time t — and decays in time at rate T + T12 
for the symmetric state and at rate F — F12 for the antisymmetric state. 

The quantity C} , defined by (|15|l . associated with C2 is then given by 

C^it) = \/[Pssit) - Paa{t)f - [psait] - Pas{t)f + 2^ p,,{t)pgg{t) , (34) 

and the negativity A/", given by 114|) , evolves in time according to the formula 

J\f{t) = max{0,A/'2(t)} , (35) 

where 

N2 {t) = ^C2{t)Ct{l^ + [Pgg{t)+Pee{t)f ~ [Pgg (t) + Pee (t)] • (36) 

Equations H33|l and 1)36(1 are exact analytical formulas describing the time evolution of entanglement 
created in the system of two identical atoms via the process of spontaneous emission. 

Let us now consider two special cases of the initial conditions: (i) initially only one atom excited, 
P44(0) = 1, (ii) both atoms initially excited, pee ~ 1- 

For case (i), we have pss{0) — Paa(O) = Pas{0) — Psa{0) = 1/2, Pee{0) — 0, and equation (p!!^ takes 
the form 

C2{t) = iv'[e-(r+ri2)t_e-(r-ri2)t]2_^e-2r*sin2(2f]i2i) . (37) 

From 1(37(1 it is seen that €2(0) = 0, there is no entanglement a,t t — 0, as it should be since the initial 
state is a product state. However, for t > 0, C2{t) becomes positive, which means that the two atoms 
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become entangled, and the degree of entanglement measured by the concurrence is given by H37|l . For 
long times, all terms in (|37|l decay to zero, and the concurrence goes to zero. One more interesting 
feature of the evolution is seen from (|37|l , it is the oscillatory behaviour of the concurrence which can be 
observed at times shorter than (2r)^^, when the oscillatory term contributes significantly to the evolution. 
The oscillations are with the frequency 2fli2, which is equal to the separation of the symmetric and 
antisymmetric states, and the oscillations become visible when the dipole-dipole interaction is sufficiently 
strong, i.e. for ^ll2 ^ F. For times longer than (2F)~^ the only term that survives is the term that 
decays with the rate F — F12, which comes from the evolution of the slowly decaying antisymmetric state, 




Figure 2. Time evolution of the concurrence C{t) (solid line), paa{i) — pss(t) (dashed line), and 
paa(t) + pss{t) (dashed-dotted line) for initially one atom excited (^44(0) = 1) with p, ± ri2, and 
ri2 = A/12 (ri2 = 0.95r, f7i2 = 9.30 T). 



and the concurrence becomes equal to the population of this state. We have numerically studied this 
behaviour, but for the interatomic distances not so short as to reveal the oscillations in the concurrence. 

In Fig.[21wc present the oscillatory behaviour of the concurrence (|37|) for the case of initially one atom 
excited and the interatomic distance ri2 = A/12, which gives the values F12 = 0.95 F and flu — 9.30 F. 
The envelops of the oscillations are given by paa [t] — Pss [t) for the lower envelope and paa (t) + Pss {t) for 
the upper envelope. The value of concurrence at the first maximum is 0.86, which is quite remarkable. 
After the time t ~ (2F)~^ when the symmetric state is practically depopulated, the concurrence becomes 
equal to the population of the antisymmetric state. 

The negativity for this case, from (|34|I - H36|I . takes the form 

A/'2 it) - { [e-(r+ri.)* _ e-(r-r,.)t] ^ + ^-^r* sm'(2n^2t)} + pl^{t) - p,, (t) , (38) 
where 

Pgg(t) = l-i[e-(r+ri2)t + e-(r-ri2)t] . (39) 

In Fig.|31we compare the time evolution of the two measures of entanglement: concurrence and negativity 
for the same values of the parameters as in Fig. [3 Generally, the negativity takes smaller values than the 
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Figure 3. Comparison of concurrence (solid) and negativity (dashed) for the same parameters as in 

Fig. El 



concurrence, except for the initial value which is zero for both of them, and the value for t 
is also zero. 

For case (ii), we have Pee(O) = 1 and the concurrence H33|l takes the form 



C2(t) = 

with 



r + r 



12 



-(r+ri2)t 



,-2rt 



r-r 



12 



r-r 



12 



r + ri2 



,-(r-ri2)t 



-2rt 



2e 



-rt 



oo which 



(40) 



Paait) = 1 - 



r + ri2 / -(r 



+ri2)t 



-2rt 



r-ri2 



r-r 



12 



r + ri 



-(r-ri2)t 



-2rt 



-2rt 



(41) 



Again, for t = the concurrence is zero, but now it is not easy to see if C2{t) can be positive, and 
numerical evaluation is needed to check the positivity. What is clear from (|4U|) . however, it is the fact 
that there are no oscillations in this case. One can also expect that for times for which the populations 
of the excited state and the symmetric state, which decay much faster than the antisymmetric state, are 
already close to zero, it is still some population in the antisymmetric state and C2{t) becomes positive. 
Numerical evaluation of l|40(l confirm that it is really true. Corresponding formula for the negativity 
can be obtained from and but this are just simple substitutions, so we do not write it 

explicitly. In Fig. 21 we plot the time evolution of the concurrence C{t), the negativity Af{t), and the 
population Paa{t) of the antisymmetric state for the initial state of both atoms excited (pee(O) = 1) and 
the interatomic distance ri2 = A/12, which gives the collective damping ri2 = 0.95 r and the dipole- 
dipole interaction frequency ^li2 = 9.30 r. As expected, there is no entanglement before the populations 
of the exited state and the symmetric state depopulate, but some entanglement appears for longer times, 
and the concurrence again becomes equal to the population of the antisymmetric state. The values of 
the negativity in this case are much smaller than the values of the concurrence, which itself is very small. 
Exciting two atoms initially is thus very ineffective in producing entanglement. 

For nonidentical atoms, although the analytical solution is possible, the formulas are rather lengthy 
and we will not adduce them here. Instead, we plot in Fig. an example of the evolution for the 
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concurrence for the case of atom "1" excited (p44(0) = 1) with A = {lj2 — wi)/2 = lOF and ri2 = A/12 
(ri2 = 0.95 r, fli2 = 9.30 r). This means that we have A ^ 1712. As it is evident from the equations of 
motion (|27|l . for two nonidentical atoms, i.e., for A ^ 0, there is a couphng between the populations of the 
symmetric state and the asymmetric state, which leads to a coherent transfer of population from one state 
to the other. This introduces oscillations into the populations of both states. Since the oscillations of the 
two populations, for the parameters of Fig.|Sl are opposite in phase, they add up in Paa{t) — Pss{t) and 
subtract in Paa{t) + Pssit), as clearly seen from the figure. The concurrence is oscillating, similarly to the 
situation shown in Fig.|21 between the lower bound {paa(t) — Pssit)) and the upper bound Paa{t) + Pssit), 
but this time the lower bound itself undergoes oscillations, which results in increasing the concurrence at 
the maxima. The value at the first maximum is 0.88, which is higher than the corresponding value for 
identical atoms equal to 0.86. It is thus possible to enhance the transient entanglement in the two-atom 
system when the two atoms are nonidentical. 

5. Conclusion 

In this paper we have studied entanglement created in a system of two two-level atoms via the spontaneous 
emission. We have found analytical formulas for the concurrence and the negativity, the two different 
measures of entanglement usually used to quantify the amount of entanglement. Our formulas are valid 
for a broad class of initial conditions which are represented by the block form of the density matrix. 

We have shown that for short times, when initially only one atom is excited, the amount of 
entanglement exhibits oscillatory behaviour until the time at which the population of the symmetric 
states becomes zero. For long times the concurrence becomes equal to the population Paait) of the 
slowly decaying antisymmetric state. For both atoms initially excited there are no oscillations, and 
the entanglement appears only for long times, when only the antisymmetric state contributes to the 
entanglement. 
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